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Abstract. In this series of lectures it is illustrated how one can study the strong dynamics of 
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of the system and the operators relevant for the reaction. Emphasis is put on the electromagnetic 
interaction with few-nucleon systems. The Lorentz integral transform method to calculate the 
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experimental results are shown. The dependence of the response functions on the nuclear interaction 
and in particular on three-body forces is emphasized. 
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Emphasis is put on the electromagnetic interaction. 



INTRODUCTION 

Electroweak (e.w.) probes are essential to study the structure and the dynamics of nuclei. 
They help to shed light on the most fundamental nuclear physics issues. In particular they 
allow 

• to assess the relevant degrees of freedom (d.o.f.) describing a nucleus. Traditionally 
these d.o.f. are supposed to be baryons and mesons. They are called "effective" 
d.o.f. to distinguish them from the fundamental d.o.f. of the strong interaction, i.e. 
quarks and gluons. Only in certain conditions of energy and momentum transferred 
by the e.w. probe to the nucleus they become relevant to interpret experimental 
results. These traditional effective d.o.f. of nuclear physics can also be divided 
in explicit and implicit ones. The former are those appearing in the Hamiltonian 
explicitly, namely the protons and the neutrons, the latter are hidden in the potential. 
A classical example of the latter d.o.f. is the pion implicit in the One-Pion-Exchange 
Potential (OPEP). 

• to assess the potential model. Nowadays various realistic nucleon nucleon (NN) 
potentials are available, that reproduce thousands of NN scattering data with very 
high precision. While they are all equivalent in describing the strong NN reaction, 
they are not, in principle, in describing a nuclear system undergoing an e.w. reac- 



tion. Therefore, by means of an e.w. probe one hopes to be able to discriminate 
among them, and have a better information about their origin. Moreover, since the 
nuclear potential has an effective nature operating between composite systems, it is 
in principle a many-body operator. An aspect of nuclear dynamics that has attracted 
a lot of interest in the last years is the importance of multi-nucleon forces and in 
particular of the three-nucleon force (3NF). Electroweak reactions can shed some 
light on their importance and origin. 
• to help understanding, by comparing theory and experiment, the microscopic ori- 
gin of typical many-body phenomenologies, like for example collectivity, clusteri- 
zations or typical single particle (mean field) behaviors. 

The e.w. probe is mediated by the photon y or by the W ± ,Zq gauge bosons. The 
coupling constants are small: the electromagnetic one a is 1/137, and the weak one 
is even suppressed by a factor 10~ 5 , due to the large mass of the weak gauge bosons. 
Therefore in calculating e.w. cross sections with hadrons it is perfectly legitimate to use 
the Born approximation, i.e. the one vector-boson approximation. This has the practical 
consequence to allow to separate the known information about the e.w. interaction from 
the unknown strong one. 

The plan of these notes is the following. 
In the first section I will outline the main steps which bring to the derivation of the 
e.w. cross section with an hadron. Extensive derivations can be found in a series of 
both classical and more modern books and articles (see e.g. iflH- flOll ). I will deal with 
the electron scattering cross section as an example, suggesting how the procedure is 
modified in other lepton scattering cases. The photoabsorption cross section will be 
viewed as a particular case of the electron scattering cross section. 
The second section will be an intermezzo to recall the scope of e.w. studies and to point 
out the interesting parts of the whole formalism. 

The third section will be devoted to one of the main ingredients of the cross section i.e. 
the four-current operator and its connections to the potential. 

In the fourth section a few considerations about the importance of ab initio approaches 
and of the study of few-body systems will be made, giving a short overview of the 
theoretical problems and how they are treated in the literature. 

The fifth section will be devoted to review the Lorentz integral transform (LIT) 
method 111 in . describing also its practical implementation. 

In the sixth section an overview of interesting applications of the LIT method to electro- 
magnetic cross sections will be given, concentrating in particular on what can be learned 
about the multi-nucleon nature of the nuclear force. 

Finally it will be concluded summarizing the main messages contained in these notes. 




Figure 1. (a) Electron scattering in the one-photon exchange approximation, (b) photoabsorption. 



OUTLINE OF THE DERIVATION OF ELECTROWEAK CROSS 

SECTIONS 

General considerations 

In this section I describe how it is possible to get information on the dynamics of a 
nucleus (and of an hadron system in general) letting it interacts with an e.w. field and 
measuring the relative cross section. What I want to put in evidence in the theoretical 
derivation of this cross section are the matrix elements < fh\jh\h >• Here the initial 
and final state of the hadron are indicated by > and \fh > respectively and jh 
is the four-current density operator (excitation operator) which is responsible for the 
internal excitation of the target, via the interaction with the external field. As was already 
mentioned in the introduction I specialize to the case of electron scattering, suggesting 
how the procedure is modified in other lepton-scattering cases. The photoabsorption 
cross section will be viewed as a particular case of the electron scattering cross section. 
Let's repeat that we work within the one photon exchange approximation, which is 
depicted in figure [TJ In figure [Qa) e and e' represent the electron before and after the 
scattering, k and k' are its initial and final four- momenta. Therefore one has: 

k = (E,k) , k' = (E',V) , E 2 - \k\ 2 = E' 2 - \k'\ 2 = m] . (1) 

Momentum and energy transferred by the electron to the nucleus are indicated by q and 
CO, respectively: 

q = k-k', (0 = E-E', CO 2 - \q\ 2 = q 2 = -Q 2 . (2) 

In figure [TJa) y* indicates the so called virtual photon. Its virtuality is given by the 
fact that q 2 (i.e. the square of the mass of this exchanged photon) is not zero, and even 
negative. This means that Q 2 is positive as it is demonsrarted in the following. 

Q 2 =\q\ 2 -(D 2 = (\k\ - \k'\) 2 - {E-E'f 

= \k\ 2 + \k f \ 2 - 2k -V -E 2 -E' 2 + 2EE' 

= \k\ 2 + \k ; \ 2 - 2 \k\ \k ; \ cos 6 - \k\ 2 - \k ; \ 2 + 2 \k\ \kf\ 

= 2\k\\k'\(l -cos0) 

= 4\k\ |^| sin 2 1 >0. 



Similar pictures as in figureQJa) can be drawn for other lepton scattering processes. One 
has simply to replace the electron with the lepton of interest and the photon with the 
appropriate gauge boson. In the limit Q 2 = the photon becomes real. Therefore the 
picture describing photoabsorption becomes figure [TJb). 

The advantage of studying an electron scattering process with respect to photoabsorp- 
tion is the fact that one can vary energy and momentum transfer independently. There- 
fore one may explore the cross section in the co\q\ plane. This aspect is strictly connected 
with the possibility of investigating the dynamics of the nucleus not only at different ex- 
citations but also at different ranges, being even able to realize the composite internal 
structure of nucleons. This happens at high momentum transfer when the wavelength 
associated to \q\ is comparable or smaller than the nucleon size. The only restriction is 
the so called space-like condition Q 2 > 0. Furthermore, in contrast to hadron probes, the 
virtual photon, as the real one, has a much larger mean free path, so that it can explore 
the whole target volume. In fact hadron probes tend to interact only on the surface. 
The large flexibility of the electrons with respect to hadrons and real photons reflects 
on the structure of the cross section. As we will see this contains longitudinal and trans- 
verse components, which allow one to have different detailed informations about nuclear 
dynamics. 

Our starting point is the differential cross section defined as follows: 

\JZ if \ 2 - 
d ° = 7^£ dk'dP f p f , (3) 

" 1 \ Ji \P target 

where |^y| 2 is the probability that the system goes from the initial to the final state 
(lowest order S-matrix) i.e. |^/| 2 = h^ss 1 | < F\H int \I > \ 2 with |/ > and \F > the 
initial and final states of the whole system (incident electron + target); \f i \ is the incident 
electron flux and the sum and average over the spin states of the initial and scattered 
electron s and s' mean that we restrict to unpolarized electrons. The phase space density 
pf is connected to the differential momenta of the fragments: 

• if the hadron remains intact pf = 1 

• if the hadron breaks into 2 fragments =>- dP\ d?2 = dPfdp i.e. pf = dp 

• if the hadron breaks into 3 fragments =>- dP\dP2dP-$ = dPjdp\dp2 etc., where 
Pf is the final momentum of the hadron center of mass (cm.), while the pt indicate 
the relative momenta of the fragments. 

In the next subsection we will elaborate on the matrix element < F\Hi nt \I >, however, 
as already mentioned, what is important for us is the part which regards the target 
< fh\jh\h >> because the internal dynamics of the system is contained in that. Therefore 
in following the evolution of the formulas one has to keep the focus on that matrix 
element. 

It is interesting to note that one can perform several different kinds of experiments 
with the electromagnetic probe. One can make elastic scattering experiments where 
the nucleus is not excited and the energy transferred by the lepton is found exclusively 
in the target recoiling energy (of course this is never the case for photons). When this 
happens one has \i^ >=\fh>. Therefore in such experiments the focus is on bound state 



properties of the nucleus. If the energy transfer serves in addition to excite the nucleus 
one speaks of inelastic scattering. In this case the information contained in the cross 
section is not only on the bound state, but also on the excited states. Depending on the 
energy transfer (and of the nucleus) they can be discrete states or continuum (break-up) 
states. One can also study absorption or emission of photons. 

If the energies and momenta do not exceed a few hundred of MeV, and the nuclei 
contain a not too large number of protons a good theoretical framework for calculat- 
ing the cross sections of all these processes consists of first order (one-photon-exchange 
approximation) quantum electro-dynamics (QED) and non relativistic quantum mechan- 
ics. For heavier nuclei the first order can become questionable. Since in the following I 
will concentrate mainly on light nuclei this problem has no relevance. 



The transition matrix 

Let us now discuss the ingredients of the transition matrix =< F\H int \I > whose 
square modulus is the main term in the cross section. 
The interaction Hamiltonian is 

H int = jd 4 xj e ^x)-A^(x), (4) 

where jt (x) is the four-current density of the electron and A^ 1 (x) is the electromagnetic 
four-potential created by the target, 

A*(x) = Jd 4 yD F (x-y)jf t (y) 1 (5) 

where Dp(x — y) is the Feynmann propagator. 
The incident electron current is 

fp = -e¥(x)Y f iV(x), (6) 
where is the Dirac spinor J~^Shu{k,s)e x . Therefore one has: 



EV 

<F\H int \I> = - em u{k\ s') y M u{k, s) < e ik> - x \ f d 4 yD F {x - y)\ e - lbx > 
VEE'V J 



<f h \e' P ^f h e-' P '-y\i h >j^- 4 . (7) 



(2k) 

Performing the integrals in d 4 y and d 4 x one gets 



<F\H int \I> 



'£^-u(k' \s')y^u{k,s) [ ^8^(k' -k+q)8^(P f -P i -q)jr (8) 
EE'V J Q 



Notice: It is here that the matrix element of interest < fh\j^\ih > appears. For econ- 
omy of notations it has been denoted simply by JP. 



The integral in d 4 q gives 



- g2m ,.nj js„ ..,u ^ {A) (Pf-Pi + k' -k) 



< F\H int \I >= ^ = ^u{k',s')y ll u{k,s) ^ { '-j* , (9) 

s/eWv q 2 

where q = k — k' = Pf — Pi. 

The lepton and hadron tensors 

If one inserts this expression in the transition probability and use the following 
relation yj] 

[{Infd^tPf-Pi + k 1 -k)] 2 = d^iPf-Pi + k 1 -k)-V-T{2n) A , (10) 
one obtains 

1 2 = (iV + fc/ ~ fe) • ^ • r E I ^(^^ ^) rM^(^, ^) 1 2 ^ v • (id 

Defining the lepton tensor 

w MV = \u(k',s')Ynu(k,s)\ 2 (12) 

and the hadron tensor W^ v = J^*J V , and considering that the electron flux can be 
written as \J e \ = pvi = y£- the differential cross section becomes 

do= e ^t\y Wllv W^8^dk l dP f . (13) 
E'\k\q 4 y M ' 

Manipulating the lepton tensor by using the properties of the traces of Dirac matrices 
one has 

w MV = -^[k^+k'^ky-g^ik-k' -m 2 )}+ih£^ va pk a k' l3 } 1 (14) 

where h is the helicity of the longitudinally polarized electron. For unpolarized electrons 
such a contribution vanishes. Let's concentrate just on the part of the cross section which 
does not depend on the electron polarization, i.e. 

4 2 i 9 

do = y^^^[k^+k^k v -gp V (k-k' -m 2 )]jV*J v d( 4 UV dP f . (15) 
The use of charge conservation 



At this point one can make use of the continuity equation q^J^ = (charge conserva- 
tion) both for the electron and hadron currents. Remember that do contains essentially 



\j e j 1 J^\ l ■ The continuity equations imply — ^ and Jo = ^r- Therefore 



= foJ° ~ Je ■ J= { ^ Je) ^ J) - I ■ J = % ■ U^P- - f] =j e .f, (16) 

with T = J — q^-p-. This means that the use of the continuity equation allows to 

% 

replace the 4-indexes \i , v = 0, 1 , 2, 3 with the 3-indexes i, j = 1,2,3 and the differential 
cross section becomes 



da 



e 4 2 
E'\k\ q 4 



[kik'j + k'ikj - gij (k-k 1 - m 2 )] J n * J ,j 8 {4) dk' dP f . (17) 



Using the relation Q 2 = —q 2 = — (2m 2 — 2k ■ k') and performing the sum on i, j one 
gets 

e 4 2 



da 



E'\k\ q q 



2Re(k-f*k'-f) + ^f*J 



d^dk'dPf. 



(18) 



Longitudinal and transverse response functions 

At this point it is convenient to choose a particular Cartesian system where 



— *. 

and decompose J in 
with 



z = q\ y= -f 7 „ ; x = yx z 
\k x k'\ 



f=f L + f T , 



qj' L = (f-q)q 
(f-x)x+(f-y)y. 



Using the definition of J and the continuity equation 7o<7o = q-J one has 

Q 2 



J', 



CO\q 



(19) 



(20) 



(21) 



(22) 



Substituting in (|T8l) the target dynamics according to the decomposition (|20l) one obtains 
three different contributions to the cross section i.e. the longitudinal, transverse and 
mixed ones. Making use of a few kinematic relations one has 



(da)! 



E'\k\ q 4 \q 4 



\ (AEE' - Q 2 ) | p | 2 8 W dk' dP f . 



(23) 



Analogously the transverse and mixed contributions are given by 



(do) 



1 



E'\k\ q 4 



4\k\ \k'\ 



\k 




k' 


\sin 2 G 


2\ 


q 


2 



Q 2 



+ 



\k\ 


\k' 


\sin 2 


2\ 


q 


2 



(da] 



LT 



-e 4 4 Q 2 
E'\k\q 4 \q\ : 



+ ■ 

4\k\ \k'\ 
d^dk'dPf, 



(24) 



\k\ \k'\sinG (E + E')Re(p*f x )8 {4) dk' dP f . (25) 



The transverse contribution consists of two separated terms which are indicated with 
the labels T and TT respectively. So at the end one can summarize the cross section as 



(do) = (V L R L + V T R T + VttRtt + V lt Rlt) <5 (4) dkf dP j 



f- 



(26) 



where the V's indicate the coefficients of kinematic nature and the R's the structure 
functions containing the matrix elements characterizing the target dynamics: 



(A 1 (4EE'-Q 2 ): 



E'\k\q A \q A 



V T 



1 



E'\k\ q* 



2\k\ \k'\ 



|X||^|sin 2 Q 



+ ■ 



2\k\ \k'\ 



e 4 1 2|fc||fc / |sin 2 e 



Vi 



LT 



E'\k\ q 

-e 4 4 Q 2 
E'\k\ q 4 \q? 

Rl = 



2\q\ 



\k\ \k'\(E+E') sin 0; 



Rt — \J' X \ + \J' X \ ; 
Rtt = \J X \ — \J'x\ ' 
R LT =Re(p*J' x ). 



(27) 

(28) 

(29) 

(30) 

(31) 
(32) 
(33) 
(34) 



Equation (1261) is very general. The kinematic coefficients (|27I) -(I301) can be expressed 
in the non-relativistic or ultra-relativistic limits for the electron kinematics. In the ultra- 
relativistic limit i.e. when the mass of the electron is much smaller that its momentum 
(that is certainly the case for electrons of tens or hundreds of MeV) E = \k\ and E' — 
then one has 

,2.. Q 4 



\k'\ V L = o M - 



(35) 



where 



1 cos 2 1 __ , 



° M = 7^2 ■ 4 fl = °Ruth. cos - ; (36) 
4£ z sin 4 1 2 

E'\k\ q 4 \ \q\ 2 2\k\\k' 

| ? |V rr 4l%^; (38) 
E'\k\ <T 2 m 

\k'\\ T = ^^^$-r\k\\k'\sme (E +E') . (39) 

E'\k\<r\qV 



The cross section for inclusive, unpolarized electron scattering 



From (1261) one can obtain inclusive, seminclusive, exclusive, etc. cross sections ac- 
cording to what one is going to measure. 

For example: if one wants the inclusive cross section i.e. one reveals the outcoming 
electron and nothing else 

— = dPf( V L R L + V T R T + V tt Rtt + V lt Rlt ) 5 (4) (P f - P, + k' -k) . (40) 
dkf J 

If one wants the seminclusive cross section i.e. one reveals the outcoming electron and 
a proton in coincidence (pf = dpp) 

-4^ = [ dP f (V L R L + V T R T + V tt Rtt + V lt R L t ) 5 (4) (P f - P x + k' - k)^- , 
dk'dPp J dP p 

(41) 

where pp is the momentum of the proton in the cm. system, while Pp is measured in 
the lab. system, etc. Let us concentrate on the simplest cross section i.e. the inclusive 
one. This is certainly the simplest from an experimental point of view, since one only 
needs (besides accelerator and target) only an electron spectrometer counting the elec- 
trons with a given energy, momentum at a certain scattering angle. A seminclusive cross 
section requires an additional hadron spectrometer and the number of counts in coin- 
cidence will be certainly smaller (smaller cross section). From the theoretical point of 
view, at a first sight, the inclusive cross section does not appear to be the simplest, but on 
the contrary in some cases the most complicate. I will comment about this in the section 
where the LIT method is described. 

A simplification of the expression for the cross section arises if one considers unpolar- 
ized targets. In fact in this case one has to sum and average on the target spin projections 
in the final and initial state, respectively. One can show that this implies that the terms 



VttRtt and VltRlt vanish. Therefore one has 



do 

dk' 



do 



da e d\k\' 

k' 2 J dP f ( V L R L + V T R T ) 5 (3) {P f - Pi - q) 8 {E) - Ef - a) . 



(42) 



where Eh , E\ represent the energies of the nucleus in the initial and final state. One 
has to stress here that these energies are not only the internal energies, but include the 
recoil energy. At this point one substitutes the matrix elements of interest in Rl and Rj 
in ([3T])-([32]) and performs the integral in dPf, obtaining 



do 



dQ. P dco 



+ 



°m ^lL\(Jn\p\ih)\ 2 8{E h f -E^-(o) 

+ 1 1 (fn I ?t I h ) \ 8 (E h f - Ef - co) 

2 9 2 J n 



(43) 



where the 5-function 8 (E h f - E\ - co) = 8 {Sf + E rec -Si - co) with Sf and S t in- 
dicating the internal energies and E rec the recoil energy acquired by the nucleus (non 
relativistically E rec = \q\ 2 /2M\ with Ma representing the mass of the nucleus). 

Notice: The integral in dPf has an important consequence: the operators p and Jj are 
function of \q\ and are expressed in the cm. system of the nucleus. 
The following remarks are in order here. 

• In order to get separate experimental information on Ri and Rj one needs to per- 
form the so called Rosenbluth separation which consists essentially in performing 
two (or more) measurements of the cross section at fixed co and \q\ and different 
scattering angles. Representing these data in a XY plot where X = Vt and Y is the 
cross section, one obtains a straight line whose properties are connected to Ri and 
Rf- 

• One can generalize the derivation of the cross section obtained above to take into 
account also the parity violating contributions. This corresponds to add to the graph 
in figure [TJ a) a similar one with the Zq boson replacing the virtual photon. 

• One can also generalize the procedure explained above to neutrino reactions (v, v') 
(V,^i ± ). In the former case the virtual photon is replaced by Zq. In the latter case 
by W ± . Of course in these cases the coupling constant a 2 = e 2 /he is replaced by 
the weak one. 

One can read more extensive derivations of the e.w. cross sections in 



INTERMEZZO 



At this point, it is better to remind what is the scope of our study. In general one could 
distinguish two different attitudes in studying e.w. interactions with nuclei. One is what 



I will call the "service" attitude. It means that nuclear theorists calculate e.w. cross 
sections that are important to solve problems of astrophysical relevance. For example 
in figures [2] and [3] one can see two famous nucleosynthesis cycles for the production of 
4 He. One can notice how many e.w. reactions are involved in the cycles. To know the 
different reaction cross sections is important to explain the abundances of elements in 
the Universe. Other e.w. cross sections are fundamental to explain the stellar evolution. 
Many of these cross sections cannot be measured in the laboratory since it is often not 
possible to reproduce the astrophysical conditions. Therefore nuclear theorists try to help 
to estimate those cross sections, as well as they can, using, in most cases, their experience 
with models or experimental inputs from other sources or, in rarer cases, accurate ab 
initio calculations. The other attitude, that I would describe as a more "fundamental" 
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Figure 2. The proton-proton cycles for the synthesis of 4 He 
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Figure 3. The Carbon-Nitrogen-Oxygen cycle for the synthesis of 4 He 

attitude consists in trying to understand what are the fundamental degrees of freedom 
(implicit and explicit) at the nuclear scale; what are the properties of the nuclear force 
and try to connect them to the underlying fundamental d.o.f. of QCD; what is the role 
of symmetries and their microscopic origin; what is the microscopic origin of nuclear 
phenomena. All that requires ab initio approaches. It means that one has to solve the 
many-body problem as accurately as possible, controlling the errors. Only in this case 
the comparison of a theoretical result, obtained with a certain input, and the experimental 



data gives information about how good the input is, without the uncertainties due to 
approximations in the solution of the many-body problem. 

In deriving the inclusive unpolarized electron scattering cross section we have seen 
how the interesting and fundamental nuclear physics ingredients enter the cross section. 
They are contained in \i^ ) and (fy | which are eigenfunctions of the nuclear Hamiltonian. 
Therefore they will be different for different nuclear potentials. However, there is more 
information in the cross section: in fact the nuclear charge and nuclear current operators 
appear. This constitutes an important difference between e.w. and purely hadron reac- 
tions. In the latter no information about charge and currents is present. This is relevant, 
since it gives access to what one calls underlying or implicit degrees of freedom. Let us 
explain this further in the next section. 



THE NUCLEAR FOUR-CURRENT OPERATOR 

If one thinks that the most reasonable effective d.o.f. in a nuclear Hamiltonian are 
protons and neutrons it becomes natural to describe the interaction of the electron with 
the nucleus as a sum of the contributions of all the interactions with the nucleons, 
thought as Dirac particles with anomalous magnetic moments. However, since one then 
calculates the matrix elements within non relativistic quantum mechanics one has to 
perform a non relativistic reduction of this interaction. This can be done via the so 
called Foldy-Wouthousen transformation [ IHQjl]. Essentially what one gets from such 



a transformation is an expression for the charge and the current operators at lowest 
relativistic order. The result is very intuitive. For the density operator one obtains 

z 

P{4) =e£exp(j#-r ( ). (44) 

This corresponds to the Fourier transform of the sum of Z 5-functions centered in the 
positions of the protons (as was already noticed above the positions have to be taken 
in the cm. of the nucleus). In this expression e is the proton charge. The neutron has 
no charge so it does not appear. Of course this ansatz is too drastic in that one knows 
that protons and neutrons have form factors. So what is done in practice is to replace 
e with the proton form factor G P E and sometimes even add the terms with the neutron 
contributions via G' E , even if in most cases they give negligible contributions. 

At lowest order the expression for the transverse current is also intuitive. One gets two 
different currents, one due to the motion of protons, called convection current and the 
other due to the anomalous magnetic moments /!,- of both protons and neutrons, called 
spin current; 

7 r -» ^ 

$ xr f Pi a* -a I 



However, this is not all. In fact one has to notice that p,7 and the Hamiltonian H are not 
independent. They have in fact to satisfy the continuity equation expressed by 

V-j = p]. (47) 

So H and p fix the divergence of the current, but not the curl, therefore they do not fix 

J. The convection and spin currents in (1431) and (|46l) are one body operators. It is easy to 

prove that only the commutator [7\p] actually reproduces the divergence of J, namely 
one has 

V-(f c +J s ) = -i[f,p}. (48) 

But then there must exist another current which is a two-(or many-)body operator which 
satisfies 

V-f exc = -i[V,p]. (49) 

Therefore one has that a given potential can only fix the divergence of this current. But 
what about the curl? 

If the potential is based on meson theory, one of course knows the entire current. 
This is the current of the exchanged meson. Since this meson does not appear as an 
explicit degrees of freedom in the Hamiltonian, it is an implicit d.o.f.. In this sense one 
says that the e.w. interaction gets information on implicit degrees of freedom, namely 
those intermediate ones, connecting nuclear physics to the fundamental theory of strong 
interaction, i.e. QCD in the non perturbative regime. 

Contrary to potentials based on meson theories, phenomenological potentials are not 
built knowing the underlying degrees of freedom. Therefore such potentials may be 
very accurate in reproducing nucleon-nucleon scattering data, but their reliability in 
electromagnetic interaction is in principle unknown since the exchange currents are not 
known. As already said charge conservation can give the constraint on the divergence 
of the currents, but no constraint on the curl exists. Therefore in order to calculate an 
electromagnetic cross section one must have a "model" for the underlying degrees of 
freedom in the potential and the comparison with electromagnetic data will allow to 
judge its reliability. 

AB INITIO APPROACHES AND FEW-BODY PHYSICS 

As already explained above ab initio calculations are those requiring the Hamiltonian H, 

the four-current (p , J) (provided the consistency in (l47l) ) and the kinematic conditions of 
the reaction as only inputs, treating all degrees of freedom of the many-body system 
explicitly and accurately (microscopic approach). Only in this way the comparison 
theory-experiment can be meaningful regarding the reliability of the inputs. However, 
ab initio approaches are a real challenge. At present only when the number of particles 
is relatively small one is able to get accurate solutions of the quantum mechanical many- 
body problem, without the need of approximations. They are necessary and unavoidable 
for more complex systems. From here comes the importance of few-body physics. 



Table 1. Kinetic (T), potential (V), binding energy £/, (all in MeV), 
and mean square radius of 4 He as obtained by various methods. 
From [32] 



Method 


(T) 


(V) 


E h 


( r 2)l/2 


FY 


102.39 


-128.33 


-25.94(5) 


1.485 


CRCGV 


102.25 


-128.13 


-25.89 




SVM 


102.35 


-128.27 


-25.92 


1.486 


HH 


102.44 


-128.34 


-25.90(1) 


1.483 


GFMC 


102.3(10) 


-128.25(10) 


-25.93(2) 


1.490(5) 


NCSM 


103.35 


-129.45 


-25.80(20) 




EIHH 


100.8(9) 


-126.7(9) 


-25.944(10) 


1.486 



Solving the quantum mechanical many-body problem has very different degrees of 
difficulty, depending if one deals with bound or continuum states. In general the situation 
is particularly problematic in nuclear physics, due to the complicated structure of the 
nuclear potential (a brief discussion about the nuclear potential will be found later). In 
the following the problems of bound and continuum states are discussed separately. 



Bound states 



Nowadays three- and four-nucleon bound states and binding energies can be calcu- 
lated with different methods based on any of the most modern high-precision NN-forces 
with an accuracy on the percentage level or less. For A=3 and 4 a well founded for- 
mulation, the Faddeev-Yakubovsky scheme (FY) [H], [TsL Oil, opened that avenue fol- 
lowed by alternative, equally accurate procedures: expansions in hyperspherical har- 
monics (HH) [fivl. ll8L Il911 orgaussians (CRCGV) [|20. l2l|l . stochastic variational method 
(SVM) ll22ll . and path integral techniques in form of the "Green's Function Monte 
Carlo" method (GFMC) [23 , 24 j25ll . Other very promising methods, based on the 



theory of effective interactions |26tl. have been developed: the "no-core shell model 
(NCSM) S Sill and the "effective interaction HH" (EIHH) HED] using expan- 
sions in harmonic oscillator and HH basis functions, respectively. 

An example of the degree of accuracy reached by these methods in the four-body 
system is given in table Q] I32D, where binding energy, expectation values of kinetic 
and potential energies and mean square radius of 4 He are shown as they result from 
calculations of seven different techniques based on the same potential model. By the 
way, remembering that the experimental value of 4 He is 28.3 MeV one can say that the 
most precise techniques lead to a clear cut answer: the potential used in Q32Q underbinds 
the a -particle significantly. This is no accident. It turns out that all modern realistic NN- 
forces underbind light nuclei significantly, showing the importance of three-body force 
(see below). 

While very light nuclei (A = 2 — 4) are privileged systems for the study of fundamental 
issues like properties and origin of the strong force, as is for example the existence of 
the multi-nucleon forces, the ab initio study of bound state properties of heavier nuclei 
has its own merit with respect to new many-body phenomena like, e.g., clusterization. 
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Figure 4. Upper panel: calculated spectra of A=6-8 nuclei from IB5I1 : lower panel: calculated density 
contours of the 8 Be ground state in the lab frame(left) and the intrinsic frame (right), labeled with densities 
in fm~ 3 . 



Furthermore, one can expect that in heavier systems some of the interaction phenomena 
found in light nuclei may be modified or amplified in view of the higher average nucleon 
density ("medium effects"). 

Going beyond four-body nuclei, the low-lying spectra for up to A=8 nucleons is rather 
well described [|34J], as shown in tabled Also these results clearly show that, in relation 
to the most modern NN-forces, 3N-forces are unavoidable in order to describe binding 
energies and low-lying excited states of light nuclei. Since the number of nucleon triplets 
overtakes more and more the number of nucleon pairs with increasing A, it is clear that 
3N-forces have to be included as well in any realistic description of complex nuclei. 

For A = 4, 6, 8, . . . one may already observe phenomena which are precursors of the 



Table 2. Experimental and GFMC energies (in MeV) of particle-stable or narrow-width 
nuclear states. Monte Carlo statistical errors in the last digits are shown in parentheses. AVI 8 
indicates the phenomenological NN potential [33], while IL stands for a version of the Urbana 
3NF. FromiH. 



Nucleus 


AV1 8 
AV lo 


TT 


n,xp 


Nucleus 


AV1 S 
AV 18 


TT 

1J_j 


n,xp 


3 H(i + ) 


-7.61(1) 


-8.44(1) 


-8.48 


? Li(H 


-31.1(2) 


-39.0(2) 


-38.77 


3 He(| + ) 


-6.87(1) 


-7.69(1) 


-7.72 


7 Li(| ) 


-26.4(1) 


-34.5(2) 


-34.61 


4 He (0+) 


-24.07(4) 


-28.35(2) 


-28.30 


8 He (0+) 


-21.6(2) 


-31.9(4) 


-31.41 


6 He (0+) 


-23.9(1) 


-29.3(1) 


-29.27 


8 Li (2+) 


-31.8(3) 


-42.0(3) 


-41.28 


6 He (2+) 


-21.8(1) 


-27.4(1) 


-27.47 


8 Li(l+) 


-31.6(2) 


-40.9(3) 


-40.30 


6 Li(l+) 


-26.9(1) 


-32.0(1) 


-31.99 


8 Li (3+) 


-28.9(2) 


-39.3(3) 


-39.02 


6 Li (3+) 


-23.5(1) 


-29.8(2) 


-29.80 


8 Li (4+) 


-25.5(2) 


-35.2(3) 


-34.75 


7 He (§~) 


-21.2(2) 


-29.3(3) 


-28.82 


8 Be (0+) 


-45.6(3) 


-56.5(3) 


-56.50 


7 Li (| ) 


-31.6(1) 


-39.5(2) 


-39.24 


8 Be (1+) 


-30.9(3) 


-38.8(3) 


-38.35 



above mentioned many-body phenomena. To illustrate this point, in the upper panel 
of figure |4] we show the spectrum of 8 Be as obtained in a microscopic "Variational 
Monte Carlo" - "Greens Function Monte Carlo" calculation (VMC-GFMC), based on 
one-body orbitals with four nucleons in an a-core coupled to (A-4) one-body (£ = 1) 
wave functions. 



Figure [5] shows NCSM results [36] for larger systems. Even if they are not fully 
converged they show the power of the method. The Coupled Cluster approach 11371. l38ll 
is also very promising to explore the medium mass region. 
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Figure 5. States dominated by p-shell configurations for 10 B, n B, 1Z C, and U C. From [3 
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Continuum states 

In general, if the reaction implies a state belonging to the continuum spectrum of 
H, the challenge may become enormous, since one has to deal with the many-body 
scattering problem, which may lack a viable solution already for a very small number 
of constituents in the system. The difficulty in calculating a many-body cross section 



involving continuum states can be realized, if one considers that, at a given energy, 
the wave function of the system may have many different components (channels), 
corresponding to all its partitions into fragments of various sizes. Already in a rather 
small system of four constituents, at energies beyond the so-called four-body break- 
up threshold, the two-, the three- and the four-body break-up channels contribute. In 
configuration space the task consists in finding the solution of the four-body Schrodinger 
equation with the proper boundary conditions. It is just the implementation of the 
boundary conditions for a continuum wave function which constitutes the main obstacle 
to the practical solution of the problem. In fact, the necessary matching of the wave 
function to the oscillating asymptotic behavior (sometimes even difficult to be defined 
unambiguously) is not feasible in practice. In momentum space the situation is as 
complicated. The proper extension of the Lippmann-Schwinger equation to a many- 
body system has been formulated long ago with the Faddeev-Yakubowski equations [|39i 
l40ll . However, because of the involved analytical structure of their kernels and the 
number of equations itself, to date it is impossible to solve the problem directly with 
their help, at energies above the four-fragment break-up threshold, even for a number 
of constituents as small as just four. 



INTEGRAL TRANSFORM METHODS 

Alternative approaches to the quite challenging problem of the dynamics in the contin- 
uum are provided by integral transform methods. The Lorentz Integral transform (LIT) 



method Q4JJ] is the natural extension of an original idea [42, 11] to calculate reaction 
cross sections with the help of integral transforms. This kind of approach is rather un- 
conventional. It starts from the consideration that the amount of information contained in 
the wave function is redundant with respect to the transition matrix elements needed in 
the cross section. Therefore, one can avoid the difficult task of solving the Schrodinger 
equation. Instead one can concentrate directly on the matrix elements. With the help of 
theorems based on the closure property of the Hamiltonian eigenstates, it is proved that 
these matrix elements (or some combinations of them) can be obtained by a calculation 
of an integral transform with a suitable kernel, and its subsequent inversion. The main 
point is that for some kernels the calculation of the transform requires the solution of a 
Schrodinger-like equation with a source, and that its solutions have asymptotic condi- 
tions similar to a bound state. In this sense one can say that the integral transform method 
reduces the continuum problem to a much less problematic bound-state-like problem. 
The form of the kernel in the integral transform is crucial. The reason is that in order to 
get the quantities of interest the transform needs to be inverted. Since it is normally cal- 
culated numerically it is affected by inaccuracies, and inverting an inaccurate transform 
is somewhat problematic. Actually, when the inaccuracies in the input transform tend to 
decrease, and a proper regularization is used in the course of inversion, the final result 
approaches the true one for various kernels |43|]. However, the quality of the result of 
the inversion may vary substantially according to the form of the kernels, even for inac- 
curacies of similar size in the transforms. In particular, when, for a specific kernel, the 
accuracy of the transform is insufficient, the result may be corrupted with oscillations 
superimposing the true solution. 
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Figure 6. (a): the function f(x); (b): its Laplace transform L(y) = J f(x) exp[— xy]dx 



In rf42"ll the Stieltjes kernel was proposed and its reliability was tested and discussed in 
simple model studies. Later, in a test of the method on a realistic electromagnetic cross 
section, calculated also in the conventional way for the deuteron 114411 . it was found that 
the use of the Stieltjes kernel is not satisfactory, since it leads to quite inaccurate results. 
The problem with the Stieltjes kernel can be understood if one notices that its form is 
not qualitatively different from that of the Laplace kernel. In fact it is well known that 
the problem of the inversion of a Laplace transform is extremely ill posed when the 
input is numerically noisy and incomplete |45|]. To illustrate the difficulties with the 



inversion of such kernels consider figure [61 In the upper and lower parts of the figure 
the function of interest and its Laplace transform are shown, respectively. As one can 
see the two curves do not resemble each other. The former shows a curve with a bump, 
the latter a monotonically decreasing curve. If now one thinks that the second may be 
affected by numerical errors as a result of the calculation one understands how difficult 
it may be to reconstruct the curve in figure [(Ja). The Laplace kernel in fact has spread 
the information about the curve over a large domain. 

Nevertheless the use of Laplace transforms is common in various fields of physics, 
from condensed matter to lattice QCD, and elaborated algorithms (e.g. the maximum 
entropy method [46]) are sometimes employed for its inversion. 

The problems encountered in inverting the Stieltjes as well as the Laplace transform 
has led to the conclusion that, differently from those two cases, the 'best' kernel should 
be of a finite range. Its extension should be, roughly speaking, about the range of the 
quantity to be obtained as the result of the inversion. Actually the perfect kernel would 
be a 8 — function. In this case in fact the function and its transform would coincide. How- 
ever, this clearly does not help! However, a so called representation of the 8 — function 
would have the necessary characteristics. 

At the same time, however, the transform has to be calculable in practice. In ll4lh it 
has been found that the Lorentzian kernel satisfies both requisites and the analogous test, 
as had been performed in [14411 for the Stieltjes kernel, has led to very accurate results. 



The Lorentz Integral transform (LIT) method 

In the inclusive electron scattering process discussed in these lectures the quantity of 
interest are the response functions Ri and Rj contained in the cross section. Like any 
general perturbation induced inclusive reaction (as are the e.w. reactions with hadron 
systems) they have the following structure 

r(E)=j£dY(Q\W Y )(W Y \Q!)8(E Y -E), (50) 

where |*Fy) are solutions to the dynamic equation 

(#-£ y )|¥ y )=0, (51) 

and H is the Hamiltonian of the system. The set is assumed to be complete and 
orfhonormal, 

jfdYF r )<?j\ = l. (52) 

The integration and summation here and in (1501) go over all discrete states and continuum 
spectrum states in the set. 

We suppose that the norms (Q\Q) and (Q'\Q') are finite. One has \Q) — O^q), 
\Q') = O , \ x ¥o), where |^Po) ls the initial state in a reaction (generally the ground state 
of the system undergoing the perturbation), and O, O are transition operators. In the 



electron scattering case discussed above they are both equal to the charge or to the 
current operators. Then one has 



r(E) =l£ dy(?o\d*\¥ Y )<F r \&\Wo)8(E Y -E). (53) 

Here |*Py) is a set of final states. 

As already pointed out, when the energy E and the number of particles in the system 
increase the direct calculation of the quantity r(E) becomes prohibitive. The difficulty is 
related to the fact that in these cases a great number of continuum spectrum states 
contribute to r(E) and the structure of these states is very complicated. 

The integral transform approach to overcome this difficulty can be considered as a 
generalization of the sum rule approach, since the use of the closure property of the 
Hamiltonian eigenstates plays a fundamental role. Consider for example a simple sum 
rule for the quantity (|53l) . based on the closure property (l52l) i.e. 



l£r(E)dE = (Q\Q'). (54) 

The calculation of this quantity is much easier than a direct calculation of r(E) itself, 
since it requires the knowledge of \Q) and \Q') only. This can be obtained with bound- 
state methods since it was supposed that \Q) and \Q') have finite norms. However, 
this sum rule contains only a limited information on r(E). In order to get much more 
information about it we consider instead an integral transform 

<f>(o)=j£ K(o,E)r(E)dE (55) 

with a smooth kernel K. This yields 

4>(<7) = dy(Q\V r )K(o,E r )<F r \Q?) 

= ^dy{Q\t{aM^y){^r\Q-)- ( 56 ) 
Using the closure property (1521) one obtains 

<P(o) = (Q\K(o,H)\Q')- (57) 

Therefore equation (1571) may be viewed as a generalized sum rule depending on a 
continuous parameter o. Only with a proper choice of the kernel K the right-hand side 
of (1571) may be calculated using bound-state type methods. And, as it was already said, 
even if <J>(ct) is available (1551) may not always be solved with enough accuracy to obtain 
r(E) via an inversion of the transform. 

Our choice of the kernel K(a,E) is such that both the calculation of 4>(<r) and the 



inversion of (p5l) are feasible. We choose [41] 



K(a,E) = - \-. . (58) 

v ' 1 (E-o*)(E-o) 



Notice that the energy parameters o that we consider are complex. For convenience we 
define them as 

a = E + a R + ioi , (59) 

where Eq is the ground-state energy, and Oj ^ 0, so that K(o, E) is actually a Lorentzian 
function centered on Eq + Or, having Oj as a half width 



K(a R ,ai,E) = T ^—^ ( 6 °) 



1 

(E-E -o R ) 2 + of 
Then the integral transform (1551) becomes 



L(a R , at) =^1dE /^l^.^ ■ (61) 

L f (E — Eq — Or) + Of 

Here and in the following the integral transform 4>(c) with a Lorentz kernel is denoted 
by L(or, O/). Using the definition (|53l) it is easy to show that the quantity (I6TT) may be 
represented as 

L(ct^, C7/) = (*|* / ), (62) 
where the 'LIT functions' *P and 'P' are given by 

I*) = (^-.Bo-ote-^)" 1 ^), (63) 
I*'} = {H-E -OR-io,y l &\^ ). (64) 

These functions are solutions to the inhomogeneous equations 

(H-Eo-OR-io^W) = d\W ), (65) 
(H-E -a R -iOi)\<¥') = O'|¥o). (66) 

When O' = 6, L(o) equals to Since for Oi ^ the integral in (|6TT) does exist, the 

norm of \¥) is finite. This implies that |*F) is a localize d function. Consequently, (1651) can 
be solved with bound-state type methods. Similar to the problem of calculating a bound 
state it is sufficient to impose the only condition that the solutions of (|65l) is localized. 
This means that in contrast to continuum spectrum problems, in order to construct a 
solution, it is not necessary here to reproduce a complicated large distance asymptotic 
behavior in the coordinate representation or singularity structure in the momentum 
representation. This is a very substantial simplification. 

Obviously, localized solutions to (|65l) and (1661) are unique. Once L(o) is calculated 
r{E) is obtained by inversion of the integral transform with a Lorentzian kernel (|6TT) 
('Lorentz integral transform'). The inversion of the LIT is discussed in the following. 



The inversion of the transform 



As already pointed out, the main advantage of the LIT method to study reactions is 
that one avoids to solve the many-body scattering problem. One solves instead, with 



bound-state methods, equations of the form (|65l) . The knowledge of those solutions 
leads to the LIT of the function r of interest. A crucial part of the method is then the 
inversion of this integral transform. The inversion of this integral transform has to be 
made with care, since errors in the transform can generate oscillations. To illustrate this 
let us consider a well defined r(E) to which we add a high-frequency term A^r(E). 
The latter leads to an additional AL^(or,(Jj) in the transform. For any amplitude of 
the oscillation A^L decreases with increasing Q. This means that for some value of Q. 
A^L may be smaller than the size of the errors in the calculation. Therefore in this case 
A a r cannot be discriminated. By reducing the error in the calculation one can push the 
frequency of the undiscriminated AQ.r to higher and higher values. 

One of the methods that can be adopted to invert the transform is called the regu- 



larization method [43]. This has led to very safe inversion results. Alternative inversion 



methods are discussed in [47]. They can be advantageous in case of response functions 
with more complex structures. Up to now, however, such complex structures have not 
been encountered in actual LIT applications, since the various considered r{e) have nor- 
mally (i) a rather simple structure, where essentially only a single peak of r(e) has to be 
resolved, or (ii) a more complicated structure, which however can be subdivided into a 
sum of simply structured responses, where the various LITs can be inverted separately. 
The present 'standard' LIT inversion method consists in the following ansatz for the 
response function 

Nfna x 

r(e') = £ c n Xn(e',(Xi), (67) 

n=l 

with e' = e — e t h, where e t h is the threshold energy for the break-up into the continuum. 
The Xn are given functions with nonlinear parameters a,. A basis set frequently used for 
LIT inversions is 

X n (e,<Xi) = e ai exp(-—). (68) 

n 

In addition also possible information on narrow levels could be incorporated easily into 
the set Xn- Substituting such an expansion into the right hand side of (|6TI) (here too the 
Oi dependence is omitted) one obtains 



N m 

I 

n=\ 



L (°r) = E c nXn{GR, Of) , (69) 

where 



Oft) = f °de' f Xn{e ffi 2 . (70) 

For given Of, the linear parameters c n are determined from a least-square best fit of L{or) 
of equation (l69l) to the calculated L(or) of equation (l62l) for a number of Or points much 
larger than N max . 

For every value of N max the overall best fit is selected and then the procedure is 
repeated ioxN' max = N max + 1 till a stability of the inverted response is obtained and taken 
as inversion result. A further increase of N max will eventually reach a point, where the 
inversion becomes unstable leading typically to random oscillations. The reason is that 



L(or) of equation (|69l ) is not determined precisely enough so that a randomly oscillating 
r(e) leads to a better fit than the true response. If the accuracy in the determination 
of L{gr) from the dynamic equation is increased then one may include more basis 
functions in the expansion (l69l) . 

The LIT method has to be understood as an approach with a controlled resolution. 
If one expects that r(E) has structures of width y then the LIT resolution parameter 
O] should be similar in size. Then it is sufficient to determine the corresponding LIT 
with a moderately high precision, and the inversion should lead to reliable results for 
r(E) if in fact no structures with a width smaller than Y are present. If, however, there 
is a reason to believe that r{E) exhibits such smaller structures one should reduce O/ 
accordingly and perform again a calculation of L with the same relative precision as 
before. Such a calculation is of course more expensive than the previous one with 
larger O/, but in principle one can reduce the LIT resolution parameter O/ more and 
more. The advantage of the LIT approach as compared with a conventional approach is 
evident. In the LIT case one makes the calculation with the proper resolution, while in 
a conventional calculation an infinite resolution (corresponding to 07 = 0) is requested, 
which often makes such a calculation not feasible. 

There are several tests of the reliability of the inversion. First of all, performing 
the calculation at different Oj one has to obtain the same stable result. If O/ is too 
small the solution tends very slowly to zero, therefore for 0/ > o" lin one may have 
numerical problems, turning into large errors for the LIT. As already said above, this 
will show up in unphysical oscillations in r(E) . This means that the stable result obtained 
with Oj < o" nn is the correct one, at that resolution. Another test is the control of the 
moments, in fact the moments of r(E) can be calculating both integrating r(E) or by 
expression (1571) (with K(o,H) = H a and o integer), which needs only the bound state. 



Practical calculation of the LIT 

We have seen how the LIT method reformulates a scattering problem as a 
Schrodinger-like equation with source terms which depend on the kind of reaction 
under consideration. These equations, which we will call the 'LIT equations' are essen- 
tially the same for any reaction, differing by the source term i.e. by their right hand side 
(see (1651) . (1661)). In all cases the asymptotic boundary conditions are bound-state-like. 
Consequently the solutions of these equations can be found with similar methods as for 
the bound-state wave functions. 

As we have seen above bound state solutions of the Schrodinger equation can be found 
in different ways searching for a direct numerical solution of the differential equations 
in coordinate space (or of the integro-differential equations in momentum space), or 
alternatively using expansions on some basis set of localized functions (bound-state-like 
boundary condition). These expansion methods become more and more advantageous 
with respect to the former one, with increasing particle number. For this reason, and 
because it has been used for the vast majority of the LIT applications I discuss it in the 
following. The truncation of the basis set converts the bound-state Schrodinger equation 
into a matrix eigenvalue problem and the LIT equations into a set of linear equations. 



These equations can be solved with various iteration methods and also Gauss type non- 
iteration ones. Such strategies have the drawback that one should solve these equations 
many times, as many as the number of Or values that one needs for a proper inversion 
of the transform. There are, however, two better strategies for calculating L(o). The 
first strategy, which is called the eigenvalue method involves the full diagonalization of 
the Hamiltonian matrix and expresses the LIT through its eigenvalues. This method is 
instructive from a theoretical point of view. 

Regardless of the reaction under consideration and the process that one wants to study 
the LIT method requires the calculation of the overlap (see (l62l)-(l64l) N ) 

L(g R ,Gj) = (*|* / ) 

= (Ql jS F 1 I • \ f LJ TP W 

(H-E -Or + lOi) {H — Eq — Gr — IGi) 

where \ Q') and \ Q) contain the information about the kind of reaction one is considering. 
Seeking and |*F') as expansions over N localized basis states, it is convenient to 
choose as basis states N linear combinations of states that diagonalize the Hamiltonian 
matrix. Let's denote these combinations | <Py) and the eigenvalues £y . The index N is to 
remind that they both depend on N. If the continuum starts at E = E t / 7 then at sufficiently 
high N the states |<%) having £y < E t % will represent approximately the bound states. 
The other states will gradually fill in the continuum as N increases. The expansions of 
our localized LIT functions read as 



N /, n N 



l*)=E- 
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<Pv>, (72) 



*')=L ^N ■ l<Pv)- (73) 
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Substituting (1721) and (1731) into (1711) yields the following expression for the LIT, 



L(g r , Gi) = ^ jjf—B n V + > < 74) 



and for 12') = 10 



L(Gr, Oj) = L T^v— p n V 7 ^2 • ( 75 ) 

v (£ v — Eq — Or) J rO l 

From (1741) and (f73T) it is clear that L(or,G[) is a sum of Lorentzians. The inversion of 
the LIT contributions from the states with £y < E t j t gives the discrete part of a response 
function, whereas the inversion of the rest gives its continuum part. The spacing between 
the corresponding eigenvalues with £y > E t k depends on N and in a given energy region 
the density of these eigenvalues increases with N. (Since the extension of the basis states 
grows with N, this resembles the increase of the density of states in a box, when its size 
increases.) 



The second strategy to obtain the LIT utilizes the Lanczos algorithm [48, 49] express 



ing the LIT as a continuous fraction. It turns out that for obtaining an accurate LIT only 



a relatively small number of Lanczos steps are needed. As the number of particles in the 
system under consideration increases the number of basis states grows up very rapidly 
and the Lanczos approach seems at present the only viable method to calculate the LIT. 

The motivation to use the Lanczos approach comes from the observation that from the 
computational point of view the calculation of the LIT is much more complicated and 
demanding than finding the ground state wave function of an A-particle system. In fact, 
to obtain the ground-state wave function one needs only to find the lowest eigenvector 
of the Hamiltonian matrix. On the contrary, as it is clear from (1741) and (1751) . the complete 
spectra of H over a wide energy range should be known to calculate the LIT. Therefore 
it is no surprise that the computational time and the memory needed to calculate L(o) 
have been the limiting factors in extending the LIT method for systems with more than 
four particles. It turns out that these obstacles can be overcome if the LIT method is 
reformulated using the Lanczos algorithm. Following reference | f30h . in this section it is 
shown how this can be done. 

To this end it is assumed that the source states \Q) and \Q') are real and rewrite the 
LIT in the following form 

L(<y ) = _iim J ( Q \ 1 \q>) \ . (76) 

o> L a R + iO] + E -H J 

A similar relation connects the response function r(e) to the Green's function 

r ( e ) = -Ilm( lim G(e + /7] +E ) ) ; G(z) = (<2|— ^|<2'> , (77) 

provided that z = e + ir] is replaced by Or + iOj. This is not surprising since the prop- 
erly normalized Lorentzian kernel is one of the representations of the 5-function and 



0//kL(or) — > r(oR) for <j[ — > 0. In condensed matter calculations 115 lLl52l l53h the Lanc- 
zos algorithm has been applied to the calculation of the Green function with a small 
value of 7], and its imaginary part has been interpreted as r(e) directly. This can be done 
if the spectrum is discrete (or discretized) and r\ is sufficiently small. In our case we 
have a genuine continuum problem and we want to avoid any discretization, therefore 
we calculate L(or) in the same way, i.e. with finite Oj using the Lanczos algorithm, but 
then we anti-transform L(or) in order to obtain r(e). 

APPLICATION OF THE LIT METHOD TO ELECTROWEAK 

RESPONSE FUNCTIONS 

In this section I present some selected results obtained with the LIT method for both 
the electron scattering response functions and the photoabsorption cross section. In the 
various cases the LIT equation (|65l ) has been solved with different ab initio bound state 



methods. For A> 4 the most efficient one has turned out to be the EIHH method ll30Ll31f| . 
which, till now, has allowed to reach results up to A=7. 

The selection of reselts presented in the following is done with two scopes: the 
first is to illustrate how from the comparison between theoretical ab initio results and 
experimental data one can get information on the nuclear force. The second is to show 



how, when A increases, a microscopic calculation can give rise to macroscopic features, 
giving the possibility to study the link between the two scales (see the case A=6). 

However, before entering into that discussion it is worth showing how the LIT method 
is working in systems where the direct calculation of continuum states is possible. This 
is the case of the two-body system, i.e. the deuteron. 



Results for A=2 



The very first LIT application was the calculation of the deuteron longitudinal form 
factor Ri{q, Q)) at\q\— 440 MeV/c in (4~Hl. where the LIT has been originally proposed. 
Actually it served as a test case for the applicability of the LIT approach, since one 
can calculate np continuum state wave functions explicitly. There it was shown how 
good was the agreement between the responses obtained in the two ways. Here another 
example on the total deuteron photodisintegration is shown. 
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Figure 7. Total deuteron photoabsorption cross section up to 100 MeV (a) and 10 MeV (b): LIT result 
(solid) and from calculation with explicit np final state wave function (crosses). 

Figure |7] shows the comparison of the LIT result with that of the conventional calcula- 
tion, using the proper np scattering states. Also here one observes an excellent agreement 
between the two calculations, showing the high precision that can be obtained with the 
LIT method. 

At this point let us make a digression about the operator which is very often used in the 
photoabsorption calculations, i.e. the dipole operator D z = JjZitf . Let us briefly explain 
why one uses such an operator in the (low energy) photoabsorption cross section. 

The photoabsorption cross section can be obtained by the electron scattering one in 
the limit Q 2 = (\q\ = (0). This means that the longitudinal part of the cross section 
disappears (the photon has no charge). Therefore one remains with the part where the 



transverse current operator is present (the photon fields are transverse to q). If one makes 
a multipole expansion of the current operator and works out the cumbersome formulas, 

at some point one realizes that there are terms that contain the V ■ J. This is nice because, 
at least for these terms, it is possible to use the continuity equation again, avoiding 
the knowledge of the current. So one can connect these terms to the charge density 
multipoles. The interesting fact is that the remaining terms are of higher order in \q\ 
and therefore can be neglected if \q\, (i.e. CO) is small enough. However, if CO is small 
enough, one can also neglect all the multipoles higher than the first, i.e. one remains 
with the dipole. This is what is called Siegert Theorem. The great importance of this 

theorem is the fact that one can avoid to know J ! The information about it is intrinsically 
contained in the use of the dipole operator. For our scope this has a very important 
consequence: if we use a phenomenological potential in the calculation, the comparison 
theory-experiment will tell us whether the unknown implicit degrees of freedom, implied 
by it, are the right ones. 



Results for A=3 

For the 3-body problem in the continuum one can solve the Faddeev equations and 
calculate continuum states both in the two- and three -body break up channels. Therefore 
also in this case one can benchmark traditional results with those obtained by the LIT. 



This has been done in [|54l |55Q. Reference [54] is interesting because there the LIT 
equation has been solved translating it into the Faddeev scheme, showing once more 
that it can be solved with different techniques. The excellent result of the benchmark 



from [55] is reported in figure [8] Notice that in the figure the small uncertainty due to 



the inversion is also shown. 



Now I am going to illustrate an instructive example, taken from references [56] 



and Il57ll . of how one gets information about the subnuclear d.o.f. from the calculation 



of Rj- In [56] the Bonn potential [58], which is a potential based on meson theory was 



used, while in 115711 the potential was the phenomenological AVI 8 [33]. Therefore in the 
former case the currents are known, while in the latter a recipe is proposed in order to 
construct them so that their divergence is consistent with the potential via the continuity 
equation. However, only the comparison with data can say whether they are correct. 
In analogy with the one boson exchange potential case these currents are called meson 
exchange currents and indicated by MEC. Figure [9] illustrates the importance of these 
MEC and the quality of the comparison with existing data. 

In case of 3 He one has a rather good agreement of theoretical and experimental 
transverse response functions for the two higher \q\ values. The MEC contribution is 
essential for reaching this agreement. At \q\ = 0.88 fm _1 , however, the theoretical Rj 
underestimates the data below 11 MeV. In the triton case the situation looks worse. 
Already for the two higher \q\ values one finds a slight underestimation of the data, in 
addition the discrepancy becomes even larger at the lowest \q\. One can conclude that 
the present agreement between theory and experiment is not bad, but certainly not very 
good. The interesting point is that it seems that a different nuclear force does not improve 
the situation, since the results at \q\ = 0.88 fm _1 with the BonnA potential from Ii56ll is 
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Figure 8. Comparison of the Faddeev and LIT results for the total 3 H photoabsorption cross section. 
The dots are the Faddeev results and the two curves represent the bounds for the inversion of the LIT. 
From 01 

almost identical to the AVI 8 result (the 3 H case was not considered in Additional 
currents involving the A resonance, up to now only partially considered in the literature 
for the threshold kinematics, could probably lead to a small improvement. 



Results for A=4 

Among light nuclei the a-particle ( 4 He) is of particular importance because it has 
some typical features of heavier systems (e.g. binding energy per nucleon), which make 
it a precious link between the classical few-body systems, i.e. deuteron, triton and 3 He, 
and more complex nuclei. For example in 4 He one can study the possible emergence of 
collective phenomena typical of complex nuclei like the giant dipole resonance (GDR). 
This is the famous bump which is seen in the total photoabsorption cross section of 
all nuclei and that has been interpreted as a collective oscillating motion of the protons 
against the neutrons. Furthermore 4 He is the ideal testing ground not only for the NN 
potential but also for the multi-nucleon forces and in particular of the three-body force 
(3BF). 
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Figure 9. Comparison of experimental and theoretical results for the threshold 3 He (upper panels) and 
3 H (lower panels) transverse response function RT as function of internal excitation energy !int at three 
momentum transfers q. Theoretical results with different current operators: relativistic one-body current 
(dashed), relativistic one-body current + MEC (solid), non-relativistic one-body current + MEC (dotted). 
Experimental data from ll59ll . 



Digression on the three-body force 

As was already remarked above, the nuclear potential has clearly an effective nature, 
therefore it is in principle a many-body operator. Yet the debate has concentrated for sev- 
eral decades only on the two-body part. Such a debate has taken place essentially among 
three different points of view: those based on meson theory, purely phenomenological 
ones and more recently the point of view of effective field theory. Realistic potentials 
based on all three different approaches have been constructed, relying on fits to thou- 
sands of N-N scattering data. After that precise calculations of the triton binding energy 
have demonstrated that a two-body potential is not enough to explain the experimental 
value, the same debate is taking place regarding three-body forces. However, for the de- 
termination of a realistic three-body potential or to discriminate among different models 
one needs to find A> 3 observables that are sensitive to it. One direction that has been 
followed 11341 . 13511 is to calculate accurately bound properties of nuclei of increasing A. 
In fact it has been realized that stronger and stronger discrepancies exist between the 
binding energies calculated with high precision two-body potentials and the experimen- 



tal values (see the upper part of figure H]). Another very promising direction is to study 
electromagnetic reactions in the continuum. In fact many years of electron scattering ex- 
periments have demonstrated the power of this kind of reactions, and in particular of the 
inelastic ones, because of the possibility to vary energy ft) and momentum \q\ transferred 
by the electron to the nucleus. This allows one to focus on different dynamical aspects 
at different ranges and one might find regions where the searched three-body effects are 
sizable. 

The 4 He nucleus is particularly appropriate for these studies because of the following 
considerations: i) the ratio between the number of triplets and of pairs goes like (A — 
2)/3, therefore it is double for 4 He than for 3 He; ii) theoretical results on hadron 
scattering observables involving four nucleons |[60|l as well as 4 He-N phase shifst [|6lll 
seem to imply that three-body effects are rather large. 



Digression on the experimental situation 

In order to appreciate the importance of the theoretical results presented in the fol- 
lowing it is necessary to summarize the situation of the experiments regarding electro- 
magnetic reactions on 4 He. 

Let us start from the real photon case, which has has a longstanding history. First 
experiments on the (y,n) reaction were performed about fifty years ago. In the follow- 
ing twentyfive years various experimental groups carried out measurements for both the 
(y,/?)and (y, n) reaction channels and the inverse capture reactions Dramatically con- 
flicting results have been obtained as a result of this work. The (y,/?) data were con- 
sistent in showing a rather pronounced resonant peak close to the three-body breakup 
threshold. At the same time, the (y,n) data at low energy were very spread, and mea- 
surements showed either a strongly pronounced or a rather suppressed giant dipole peak. 
In 1983 a careful and balanced review of all the available experimental data for the two 
mirror reactions was provided Ii62ll . A strongly peaked cross section at low energy for 
the (y,p) channel and a flatter shape for the (y, n) one was recommended by the authors. 
Three new experiments on the (y, p) reaction were subsequently carried out, two of them 
contradicting and one confirming the recommended cross section. Measurements of the 
ratio of the (y,p) to (y, n) cross sections in the giant resonance region were performed 
as well ll63"ll and, at variance with the cross sections recommended in rf62ll . results very 
close to unity were reported. Finally, in 1992 additional cross section data were deduced 



from a Compton scattering experiment on 4 He Il64ll . A strongly peaked cross section for 
the 4 He total photo absorption was found suggesting a (y, n) cross section considerably 



larger than the one recommended in 116211 . 

In 1996 the first theoretical calculation of the two-fragment breakup cross section 
with inclusion of FSI was performed in the energy range below the three-fragment 



breakup threshold H65Q . The semi-realistic MTI-III potential [66] was employed. The 
results of Il65ll s howed a rather suppressed giant dipole peak. The agreement with the 
(y,n) data in [67] was very good and the situation seemed to be settled. 

However, the following year a calculation of the total photoabsorption cross section 
for the 4 He up to the pion threshold was carried out 116 811 . Full FSI was taken into 



account in the whole energy range via the application of the LIT method. The four- 
nucleon dynamics was described with the same NN potential as in Il65ll . Different from 
the previous work apronounced giant dipole peak was found. These results have been 
reexamined in [|69L fvoh . A small shift of the peak position has been obtained, but the 
pronounced peak has been confirmed. 

Unfortunately, the experimental situation of the 4 He photodisintegration is not yet 
sufficiently settled. Most of the experimental work has concentrated on the two-body 
break-up channels 4 He(y, n) 3 He and 4 He(y,p) 3 H in the giant resonance region, but still 
today there is large disagreement in the peak. In fact in two recent (y,«) experiments 
117 ll LZ2D one finds differences of a factor of two. A measurement of the analog of the 
GDR in 4 He has been performed in Q73[| via the 4 He( 7 Li, 7 Be) and points to a cross 
section which is very close to the result in Q68[|. 

Considering that from the early times of atomic physics the cross section for absorp- 
tion of photons has always represented the fundamental observable to study the spec- 
trum of composite systems, the present situation is unacceptable. Therefore it is highly 
desirable that this observable attracts the interest of experimentalists. (Projects in this 
direction have been recently considered at MaxLab in Lund and HlyS at TUNL in North 
Carolina). 



Regarding the experimental situation for electron scattering, in the '80 and '90's an 
intense experimental activity has been devoted to it, in particular to inclusive electron 
scattering (denoted by (e,e')) in the so called quasi-elastic (q.e.) regime, correspond- 
ing to momentum transfers of several hundreds MeV/c and energies around the q.e. 
peak, observed at about q 2 /2m. In such conditions one can envisage that the electron 
has scattered elastically with a single nucleon of mass m. Various nuclear targets have 
been considered, from very light to heavy ones. The reason for concentrating on the q.e. 
regime has been the conviction that for such kinematics the plane wave impulse approx- 
imation (PWIA) might be a reliable framework to describe the reaction. The neglect of 
the final state interaction (FSI) has the advantage to allow a simple interpretation of the 
cross section in terms of the dynamical properties of the nucleons in the ground state. 
In particular the focus in those years were the ground state short range correlations, i.e. 
the dynamical effects on the wave function of the largely unknown repulsive short range 
part of the potential. 

It is clear, however, that it is important to clarify the reliability of the PWIA. As it will 
be seen in the following, the LIT method, which takes into account the full dynamics in 
the final state, allows to clarify this issue. 



The comparison theory-experiment for the photon case 

An important step forward on the theory side has been made by performing a cal- 
culation of the total photoabsorption cross section Oy of 4 He with a realistic nuclear 
forc e [17411 . To this end the four-body problem with the Argonne VI 8 (AVI 8) NN poten- 
tial 113311 and the Urbana IX (UIX) 3NF 112411 has been solved. The results are shown in 
figure [101 
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Figure 10. Total 4 He photoabsorption cross section: (a) a y (AV18) and C~ 19 (AV18+UIX), (b) as 
(a) but also included upper/lower bounds and various experimental data (see text), area between dotted 
lines |67l|73], dotted box JH, squares 01, and circles fvk . 



Due to the 3NF one observes a reduction of the peak height by about 6% and a shift of 
the peak position by about 1 MeV towards higher energy. Large effects of the 3NF are 
found above 50 MeV with an increase of Oy by e.g. 18%, 25%, and 35% at CO = 60, 
100, and 140 MeV, respectively. It is very interesting to compare the present 3NF 
effects to those found for oy( 3 H/ 3 He) ||75l. l55|l . Surprisingly, the peak height reduction 
is smaller for 4 He. For 3 H/ 3 He the size of the reduction is similar to the increase of E^ 
(10%), whereas for 4 He the 3NF increases E^ by 17%, but reduces the peak by only 
6% and thus cannot be interpreted as a simple binding effect. Also at higher CO there 
are important differences. The enhancement of oy( 4 He) due to the 3NF is significantly 
stronger, namely about two times larger than for the three-body case. Interestingly this 
reflects the above mentioned different ratios between triplets and pairs in three- and 
four-body systems. In figure [Tola) also ay for the semi-realistic Malfliet-Tjon (MT) NN 
potential [|68l. l69ll is illustrated. Similar to oy( 3 H/ 3 He) [p75|l one finds a rather "realistic" 
result in the giant resonance region (overestimations of the peak by about 10-15%) and 
quite a correct result for the peak position; however, at higher energy cy is strongly 
underestimated, for 4 He by a factor of three at pion threshold. In figure [TOta) data 
from [77] are also shown. They are the only measurements of oy( 4 He) in the whole 



energy range up to pion threshold. In the peak region the data agree best with the MT 
potential, while for the high-energy tail one finds the best agreement with the AVI 8 
potential. 

In figure [T0fb) our low-energy results are compared to further data, which, however 
have to be interpreted with some care, since no of them corresponds to a direct mea- 
surement of the total photoabsorption cross section: (i) in [ ]64( ] the peak cross section is 
determined from Compton scattering via dispersion relations, (ii) the dashed area corre- 
sponds to the sum of cross sections for (7, n) from rfoTll and (7, p) 3 H from IT76ll as already 
shown in [68], (iii) the data from the above mentioned recent (7, n) experiment IT7lll are 
included only up to about the three-body break-up threshold, where one can rather safely 
assume that o r — 2c(y,n) (see also |70|]), (iv) in [72] all open channels are considered. 
One sees that the various data are quite different exhibiting maximal deviations of about 
a factor of two. The theoretical oy agrees quite well with the low-energy data of [|6Vl. F761 . 
In the peak region, however, the situation is very unclear. There is a rather good agree- 
ment between the theoretical o r and the data of and (|64l] . while those of |6_7LlZ6J] are 
noticeably lower. Very large discrepancies are found in comparison to the recent data of 
Shima et al. fZ2l], while a very good agreement is found with data in [73]. 

From all this long discussion it is evident that the experimental situation is rather 
unsatisfactory and further improvement is urgently needed, if one wants to clarify the 
3NF issue. In particular it is worth to stress that due to the above mentioned Siegert's 
theorem what one is testing here are also three-body exchange currents connected to the 
3NF. 



The comparison theory-experiment for the electron case 

Inelastic electron scattering off nuclei provides complementary informations on the 
nuclear dynamics. In fact varying the momentum \q\, transferred by the electron to the 
nucleus, one can focus on different dynamical regimes. At lower momenta the collective 
behavior of nucleons is studied. As \q\ increases one probes properties of the single 
nucleon in the nuclear medium and its correlations to other nucleons from long- to 
short-range. Since the longitudinal response is not sensitive to meson exchange 
effects (only at lowest relativistic order in the Foldy-Wouthuysen transformation) the 
use of a simple one-body density operator allows to concentrate on the nuclear dynamics 
generated by the potential and one might find regions where the searched three-nucleon 
effects are sizable. 

Here results on Ri are presented, focusing on the evolution of dynamical effects 
as the momentum transfer decreases. In figures [TT] and [121 Rl at constant \q\ = 200 
and 100 MeV/c are shown. One has a large quenching effect due the 3NF, which is 
strongest at lower \q\. One should notice that such an effect is not simply correlated to 
the under-binding of the AVI 8 potential. In fact, if this was the case, the results with the 
Malfliet-Tjon potential (MT) [66], which gives a slight over-binding of 4 He, would lay 
even below those obtained with AV18+UIX. On the contrary the MT curve is situated 
between the curves with and without 3NF. 

Given the large 3NF effect at lower \q\ it is interesting to see whether there is a 





Figure 12. Longitudinal response function for \q\ = 100 MeV/c with the AV18 (dashed), AV18+UIX 
(solid) and MT (dashed-dotted) potentials. 



dependence of the results on the 3NF model itself, investigating in addition other low-q 
values. To this end the calculation has been performed using also the Tucson Melbourne 
(TM') Q79D three-nucleon force at \q\ = 50 MeV/c. While the UIX force contains a two- 
pion exchange and a short range phenomenological term, with two 3NF parameters fitted 
on the triton binding energy and on nuclear matter density (in conjunction with the AVI 8 
two-nucleon potential), the TM' force is not adjusted in this way. It includes two pion 




Figure 13. Longitudinal response function for \q\ — 100 MeV/c with the AV18 (dashed), AV18+UIX 
(solid) and MT (dashed-dotted) potentials. 



exchange terms where the coupling constants are taken from pion-nucleon scattering 
data consistently with chiral symmetry. Figure[[3] shows that the increase of 3NF effects 
with decreasing \q\ is confirmed. Moreover it becomes evident that also the difference 
between the results obtained with two 3NF models is sizable. One actually finds that 
the shift of the peak to higher energies in the case of UIX generates for Ri a difference 
up to about 10% on the left hand sides of the peaks. This is a very interesting result. 
It represents the first case of an electromagnetic observable considerably dependent on 
the choice of the 3NF. In the light of these results it would be very interesting to repeat 
the calculation with EFT two-and three-body potentials. At the same time it would be 
highly desirable to have precise measurements of Ri at low \q\ in order to discriminate 
between different nuclear force models. 

In figure [14] an overview of the results obtained for larger \q\ is given, showing also 
the comparison with existing experimental data. One sees that the 3NF results are closer 
to the data, this is particularly evident at \q\ = 300 MeV/c. However, the 3NF effect is 
generally not as large as for the lower momentum transfers. In some cases the quenching 
of the strength due to the 3NF is comparable to the size of the error bars, particularly 
for the data from Il80n . The largest discrepancies with data are found at \q\ = 500 
MeV/c. While the height of the peak is well reproduced by the result with 3NF, the 
width of the experimental peak seems to be somewhat narrower than the theoretical 
one. On the other hand one has to be aware that relativistic effects are not completely 
negligible at \q\ = 500 MeV/c. They probably play a similar role as found in the electro- 
disintegration of the three-nucleon systems (see e.g. [|83j]). In the case of \q\ — 250 
MeV/c the experimental results are not sufficiently precise to draw a conclusion. 

But the most striking message one gets from those result is the large FSI effects that 
one finds in all cases, an effect that is essential for reaching agreement with experiment. 
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Figure 14. R L (co,q) at various \q\: PWIA (dotted); full calculation with AV18 (dashed) and AV18+UIX 
(solid). Data from Bates 08OI1 (squares), Saclay llSltl (circles) and world-data set from [82] (triangles). 



The PWIA results fail particularly in the q.e. peak and at low energies. In fact while the 
FSI effect decreases with increasing \q\ in the peak region, this is not the case at lower 
energies. This is particularly a bad news if one considers that the high |g|-low CO region 
was considered the best one to focus on ground state short range correlations. 



Table 3. Temperature averaged neutral current inclusive inelastic 
cross-section per nucleon (in 10 _42 cm 2 ) as a function of neutrino tem- 
perature (in MeV). From 1 84] 
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Inelastic Neutrino Reactions 

Here I report an example of weak cross section, calculated ab initio by the LIT 
method, which serves to clarify a problem of astrophysical interest. In fact the current 
theory of core collapse supernova holds some open questions regarding the explosion 
mechanism and late stage nucleosynthesis. In particular, due to the high abundance 
of a particles in the supernova environment, the inelastic neutrino- 4 He reaction is of 
particular relevance. The characteristic temperatures of the emitted neutrinos are about 
6-10 MeV for v M;T (v M)T ), 5-8 MeV for v e , and 3 - 5 MeV for v e . 

In |[84ll a full ab-initio calculation of the inelastic neutrino- 4 He reactions has 
been considered for the channels 4 He(v x ,v x )pi, 4 He(v x ,\'^)pi, 4 He(v e ,e + ){X, and 
4 He(v e ,e~) 4 X, where x = e,\x,x and ^X stands for the final state A-nucleon system, 
with charge Z. 

Table [3] gives the temperature averaged total neutral current inelastic cross-section as 
a function of the neutrino temperature for the AV8', AVI 8, and the AV18+UIX nuclear 
Hamiltonians and for the AV18+UIX Hamiltonian adding the MEC. From the table it 
can be seen that the low-energy cross-section is rather sensitive to details of the nuclear 
force model (the effect of 3NF is about 30%). This sensitivity gradually decreases with 
growing energy. In contrast the effect of MEC is rather small, being on the percentage 
level. 

The overall accuracy of these results is of the order of 5% and mainly due to the 
strong sensitivity of the cross-section to the nuclear model. The numerical accuracy 
of the calculation is of the order of 1%. Considering this, one can say that, thanks to 
the power of the LIT method, an important step has been done in the path towards a 
more robust and reliable description of the neutrino heating of the pre-shock region in 
core-collapse supernovae, in which 4 He plays a decisive role. 



Results for A=6,7 

Increasing the number of nucleons one may hope to find the surge of typical col- 
lective effects. This is indeed what happens if one study the total photodisintegration 




Figure 15. Total photoabsorption cross sections for the six-body nuclei with AV4\ MN and MTI-III 
potentials: 6 Li (a), 6 He (b). 



cross section of the 6-body nuclei 6 Li and He. In figure \\5\ the results for the total 
photoabsorption cross section of 6 Li and 6 He 11851. 18611 with several semirealistic poten- 
tials are shown. One notes that the general structure of the cross section is similar for 
the various potential models. In particular one has always the presence of two peaks 
in 6 He, even if peak positions and peak heights are potential dependent, while there is 
one single giant dipole resonance peak in 6 Li. The double peak structure of 6 He can 
be interpreted as a response of a halo nucleus, where the low-energy peak is due to the 
halo-a core oscillation (soft dipole response) and the peak at higher energies due to the 
neutron-proton spheres oscillation (Gamow-Teller mode or hard dipole response). So 
the low-energy 6 He peak is due to the breakup of the neutron halo. The second one, at 
higher CO, corresponds to the breakup of the a core. The 6 Li cross section does not show 
such a substructure. This is probably due to the fact that the breakup in two three-body 
nuclei, 3 He + 3 H, fills the gap between the halo and the a core peaks. Note that in case 
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Figure 16. Theoretical and experimental photoabsorption cross section results (see also text). 



of 6 He a corresponding breakup in two identical nuclei, 3 H + 3 H, is not induced by the 
dipole operator. 

In figure [16] the theoretical results are shown together with available experimental 
data. For the AV4' potential [87] active also in P-wave, one finds an enhancement of 
strength in the threshold region compared to the S-wave potentials. It is evident that 
the inclusion of the P-wave interaction improves the agreement with experimental data 
considerably. This is particularly the case for "Li. In fact with the AV4' potential one 
has a rather good agreement with experimental data up to about 12 MeV. In case of 
6 He the increase of low-energy strength is not sufficient, there is still some discrepancy 
with data. Probably, in order to describe the halo structure of this nucleus in more detail 
additional potential parts are needed. In particular the spin-orbit component of the NN 
potential could play a role in the determination of the soft dipole resonance. In fact in a 
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Figure 17. Comparison of the theoretical photoabsorption cross section calculated with AV4' potential 
with experimental data from JH]. 



single particle picture of He the two halo neutrons will mainly stay in a p- state and can 
interact with one of the core nucleons via the NN LS-force. 

Anyway the experimental situation is very confused. Again, like in the 4 He case, no 
data comes from a direct measurement of the total photoabsorption cross section. Such 
an experiment is now in course at MaxLab in Lund for 6 Li and we really hope that they 
will be able to clarify the situation. 

Finally in figure [FT] the total photoabsorption cross section of 7 Li is shown. One 
readily notes that the gross properties of the data, steep rise, broad maximum and slow 
fall off, are very well reproduced quantitatively over the whole energy region by the 
theory. It is worthwhile to emphasize that this result is based on an ab initio calculation 
in which the complicated final state interaction of the 7N-system is rigorously taken into 
account by application of the LIT method. No adjustable parameters were used, the sole 
ingredient being the AV4' NN potential model. It remains to be seen whether the slight 
variation of the data near and above the maximum will also be found in an experiment 
with improved accuracy. Such an experiment is in course at MaxLab in Lund. The new 
results could clarify in particular the question whether a simple semi-realistic potential 
like the AV4' model is sufficient for an accurate theoretical description of this reaction 
or whether a more realistic nuclear force including a 3N-force is needed. 



CONCLUSIONS 



In these lectures I have tried to give an overview of what one can learn from studying 
e.w. interactions with nuclei. The essential messages are the following. 

• Since the nuclei are interacting strongly and the e.w. interaction is weak a lowest 
order perturbative description of the cross sections allows to focus on the nuclear 
dynamics. 

• Understanding nuclear dynamics means understanding how the effective degrees of 
freedom that appear in the hamiltonian explicitly (nucleons) interact. Since these 
d.o.f. are not fundamental particles, but composite systems, the interaction is in 
principle a many-body interaction. In particular to understand the role of more than 
two-body forces in nuclei is an issue which is debated at present. 

• An effective interaction between effective d.o.f. implies implicit d.o.f.. They con- 
stitute the bridge between nuclear physics and non perturbative QCD. Electroweak 
probes, differently from the hadronic ones, are sensitive to them, providing precious 
informations about this issue. 

• In order that a comparison between theory and experiment is meaningful regarding 
the conclusions one can draw, it is necessary to work within theoretical ab initio 
approaches able to get accurate results with controlled numerical uncertainties. 
Due to the difficulties in solving the quantum mechanical many-body problem, 
especially in the continuum (many particle scattering problem), few-body systems 
assume an important role. 

• The lorentz integral transform (LIT) method has represented a big step forward 
in that it has allowed to calculate ab initio e.w. cross sections with more than 
three nucleons in the continuum, reducing the continuum problem to a bound state 
problem. 

• Until now applications of the LIT method, coupled with advanced bound state 
methods, have concentrated on systems with A=2-7. Modern realistic potentials, 
however, have been used up to A=4. The present challenge on the theory side is the 
extention to larger systems, both stable and exotic (halo). 

• Confronted with the big progresses of the theory the experimental situation does not 
seem to have reached the same amount of accuracy, hindering in many cases the 
possibility to draw conclusions from the comparison theory-experiments. Therefore 
there is a clear necessity of experimental activity in this field, especially at lower 
energies and momenta. 

In conclusion I hope to have been able to give an idea of the accomplishments in the 
field of e.w. interactions with light nuclear systems, (stable or exotic) and in particular 
of its future perspectives. They could be very exciting, considering the rich amount of 
investigations which are now made possible by the great progress in the field of few- 
body physics. 
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